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CERTAIN  PROBLEMS  OF  THE  HYDRODYNAMIC. 
OP  RELAXING  MEDIA 


I.  P.  Stakhanov  and 
Ye.  V.  Stupochenko 

Moscow 


This  article  examines  certain  general  properties  of  the  equations 
of  motion  of  thermodynamically  nonequilibrium  media  and  the  transition 
to  "equilibrium"  hydrodynamics  in  the  limiting  case  of  small  relaxa¬ 
tion  times.  Certain  problems  of  the  propagation  of  small  disturbances 
in  relaxing  media  are  examined  in  a  linear  approximation. 

In  any  motion  of  a  liquid  or  gas,  excluding  certain  special 
cases  such  as  the  uniform  translational  motion  of  a  system  as  a  whole, 
disturbances  arise  in  the  statistical  equilibrium  state  of  the 
medium.  If  the  gradients  of  the  hydrodynamic  values  are  not  too 
large,  the  equilibrium  distribution  function  is  only  negligibly 
disturbed;  however,  from  the  molecular  and  kinetic  viewpoint  these 
slight  disturbances  are  the  basis  of  all  dissipative  processes  in 
moving  gases.  Introducing  terms  with  coefficients  of  viscosity, 
thermal  conductivity,  diffusion,  etc.,  Into  ordinary  hydrodynamic 
equations  is  the  means  used  for  calculating  these  small  disturbances 
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of  BtgtlBtltll  equilibrium;  in  these  calculations  the  thermodynamic 
p trimeters  Of  the  medium  are  taken  equal  to  their  values  in  a  state 
of  complete  thermodynamic  equilibrium  (local  equilibrium). 

HoWfver,  in  a  number  of  cases  this  type  calculation  of  statis¬ 
tical  nonequilibrium  is  insufficient.  The  time  t  for  establishment 
of  local  statistical  equilibrium  varies  within  broad  limits  depending 
on  the  nature  of  the  process.  It  Is  of  the  same  order  of  time  as 
the  free  path  used  for  establishing  Maxwell  velocity  distributions 
in  a  gas  and  is  somewhat  greater  for  the  case  of  rotational  motion 
of  molecules.  Such  rapid  processes  are  outside  the  range  of  applica¬ 
bility  of  hydrodynamic  equations.  Spatial  regions  with  substantiate 
disturbance  of  statistical  equilibrium  of  this  type  are  described  in 
hydrodynamics  by  mathematical  discontinuity  surfaces.  However,  for 
the  case  of  molecular  oscillations  the  relaxation  time  is  of  the 
order  of  thousands  and  ten  thousands  of  free-path  length  times.  The 
exchange  of  kinetic  energy  of  translational  motion  between  particles 
having  great  mass  —  the  electrons  and  ions  -  is  greatly  hampered. 
Finally,  considerable  time  is  required  to  establish  chemical  or  ioniz¬ 
ation  equilibrium. 

If  the  characteristic  macroscopic  reaction  time  is  comparable 
in  order  of  magnitude  to  the  relaxation  time  x  for  any  of  the  processes 
occurring  in  the  medium,  the  thermodynamic  parameters  can  differ 
substantially  from  their  equilibrium  values,  and  the  system  of  hydro- 
dynamic  equations  should  be  supplemented  with  a  kinetic  equation 
describing  the  relaxation  process.  These  conditions  are  encountered 
in  supersonic  flows,  for  example,  in  the  zone  directly  behind  the 
front  of  a  shock  wave  which  markedly  disturbs  the  state  of  thermo¬ 
dynamic  equilibrium  [lj.  The  other  well-known  example  of  motion  in 
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which  explicit  calculation  of  the  relaxation  processes  la  necessary, 
is  the  propagation  of  ultra  high-frequency  sound.  The  studies  of 
I.  G.  Kirkwood  and  L.  I.  P.  Broer  [2,  }]  are  devoted  .0  obtaining 
a  closed  system  of  hydrodynamic  and  kinetic  equations  for  different 
processes . 

We  will  consider  below  certain  genera]  proportion:  uf  the  system 
of  "relaxation  hydrodynamic"  equations,  as  well  as  the  possibility 
of  passing  to  ohe  Jimit  to  the  usual  "equilibrium"  hydrodynamic:, 
hydrodynamics  at  small  values  of  t. 

The  latter  is  of  interest  because  certain  important  properties 
of  "relaxation  hydrodynamic"  equations  when  x  approaches  zero  are 
not  transformed  into  properties  of  the  equilibrium  hydrodynamics 
equations.  For  example,  the  characteristics  of  equations  of  unsteady 
one-dimensional  motion  lacking  viscosity  and  thermal  conductivity 
described  by  L.  I.  F.  Broer  [>]  do  not  depend  on  the  value  of  time 
x  and  do  not  coincide  with  the  characteristics  of  equilibrium 
hydrodynamics . 

As  will  be  shown  below,  this  difficulty  is  removed  if  we  assume 
that  the  equations  of  equilibrium  hydrodynamics  are  an  approximate 
description  of  only  the  special  class  of  solutions  of  relaxation 
hydrodynamic  equations.  These  solutions  are  characterized  by  the 
fact  that  local  nonequilibrium  at  every  given  moment  is  fully  defined 
by  the  fields  of  hydrodynamic  values  regardless  of  the  previous 
history  of  the  medium.  Moreover,  kinetic  equations  describing  the 
change  in  time  of  the  parameters  describing  nonequilibrium  of  the 
medium  drop  out  of  relaxation  hydrodynamics. 

We  shall  also  investigate  the  propagation  of  small  disturbances 
in  relaxing  media.  We  shall  show  that  in  the  region  of  the  order  of 


-> 
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the  "relaxation  length"  l,  they  will  be  propagated  along  the  charac¬ 
teristics  of  the  equations  of  relaxation  hydrodynamics;  therefore 
near  obstacles  weak  discontinuities  arise  whose  direction  does  not 
coincide  with  the  usually  observed  Mach  lines  and  which  exponentially 
vanish  at  distances  large  in  comparison  with  l,  At  large  distances 
the  propagation  of  the  disturbances  is  determined  by  the  equilibrium 
velocity  of  sound.  However,  unlike  equilibrium  hydrodynamics  this 
disturbance  is  made  up  of  disturbances  arriving  from  different  points 
in  this  case  as  propagation  proceeds,  the  shape  of  the  disturbances 
changes  so  that  the  "fine  details"  gradually  disappear. 

1.  Equations  of  relaxation  hydrodynamics.  We  shall  assume 
that  the  nonequilibrium  states  which  occur  are  states  of  incomplete 
statistical  equilibrium  which  can  be  described  thermodynamically. 

As  the  characteristic  function  determining  the  thermodynamic  proper¬ 
ties  of  the  medium,  we  choose  the  internal  energy  of  a  unit  of  mass, 
e(p ,  s,  0,  which  can  be  considered  as  a  function  of  density  p, 
specific  entropy  _s,  and  the  additional  parameter  Z  (or  several  such 
parameters)  [4,  5],  which  describes  the  deviation  from  complete 
thermodynamic  equilibrium.  The  physical  sense  of  the  parameters  Z 
can  be  quite  different.  This  can  be  for  example,  the  temperature 
of  che  internal  degrees  of  freedom  or  the  concentration  of  the 
reacting  component.  In  describing  the  chemical  reaction  it  is 
especially  convenient  to  Introduce  the  parameter  Z  by  means  of  the 
equation  [2] : 

dim  = 

k 

Here  dm^  Is  change  of  mass  of  the  component  i_  resulting  from 
reaction  k,  and  v^  ls  the  corresponding  stoichiometric  coefficient. 
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The  actual  system,  an  a  rule,  ia  characterized  by  ueveral  nuch  param¬ 
eters  -,  becauae  acveral  reactions  may  simultaneously  lake  place  in 
the  ayuiem,  and  moreover  both  the  temperature  of  the  dllYorent 
components  and  the  temperature  of  their  internal  (Tor  example, 
fluctuating)  degrees  of  freedom  can  differ.  We  aha  1  j  limit  ouraelvea 
to  a  consideration  of  the  simplest  one-paramei  ■;  relaxation  process, 
for  which,  however,  we  can  investigate  the  characteristic  features 
of  relaxation  hydrodynamics. 

Iii  the  state  of  equilibrium,  £  „akes  the  value  £(p,  s)  defined 
by  the  equation 

to)  =  0,  ■  e5  =  [||]ci  (1.1) 

The  equations  of  continuity,  momentum,  and  energy  have  their 
usual  form 

^  J.-  div  pv  =  0 

p(^Y  (V’V) v)  = -grudp  +  v-n  (p=Pj[|](.)  (1‘2) 

Mt  +  pe)=  _f)iv{pv(f +  <J’)- ( n'V)  +  <l'j 

Here  v  is  velocity  of  the  medium,  p  is  pressure,  u>  is  specific 
enthalpy,  q  is  the  density  vector  of  the  thermal  flow,  II  is  the  tensor 
with  components  definable  from  the  relation 

Pik  =  —  Jtj k  (  1 . 3  ) 

where  is  the  stress  tensor.  It  further  follows  that  n  need  not 
coincide  with  the  tensor  of  viscous  stresses.  The  following  form 
of  the  energy  equation  also  follows  from  (1.2): 


'/>’•  (1  1  1  i)q j  tliJj  f  d  _  0 

dl  ~  P  dl  7  _  7  ~  ^  (rTl’  ~  Hi  V 


(1.^) 
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From  here  and  from  the  general  thermodynamic  cxpressioi 


*m~*'(r}J'  »+•»*  (r-@)J 

where  T  la  temperature  of  the  medium,  we  obtain 

l/t  1  1  <>Ct  I';  l!l 

ill  (■/' i'/,  1  pi  <’ik  Till 

In  the  general  caae,  the  parameter  r,  at  a  given  point  of  the 
medium  can  vary  both  duo  to  ocular  j'low  _r  and  to  uhe  vector  flow  J. 

In  the  future  the  phyuiea.  uenao  of  the  parameter  will  not  no  speci¬ 
fied;  therefore  the  phyalca i  sense  ■  <f  flows  r  and  J  may  not  no  involved. 

In  the  particular  caae,  when  £,  defines  the  composition  of  the 
gas  dissociating  according  to  the  equation  2A 7I  A a,  flow  £  will  be 
the  dissociation  and  J  will  be  the  dii'fusion  flow. 

When  defining  r  and  J  appropriately  the  equation  for  the  rate 
of  change  of  £  can  a j ways  be  written  as 

=  ,Jivj  (1.7) 


(1.5) 


(1.0) 


Substituting  (l./)  into  (l.b),  we  obtain  the  equation  of  the 
entropy  balance 


ds 


d  ( 


pTt  +  57i('7r 


e"Ji 

T 


To  = 


1  01  ov ;  . ,  “  eE 

_  f/i  f  +  n,k  oTk~  Fxi  T 


We  easily  see  that 


d(ps) 

~7T  =  ~  <llv  0.  -I-  pi'v)  +  ct, 


(1.8) 

(1.9) 


(1.10) 


The  vector  Jt,  is  the  density  vector  of  the  entropy  flux,  or  to 
be  more  exact,  that  part  of  the  density  of  the  entropy  flux  which  is 
associated  with  macroscopic  motion  of  the  medium;  the  total  density 
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of  the  entropy  flux  i3  +  p3V. 

Tliu  right-hand  aide  of  (l.B)  represents  formation  of  entropy  o, 
i.e.,  the  rate  of  change  in  the  entropy  of  a  unit  volume  due  to  the 
irreversible  processes  oc curing  In  It:  the  thermal  conductivity, 
interna]  friction,  and  processes  dcscribable  by  flow.;  r  and  J  (in  the 
particular  ease,  the  processes  can  be  chemical  reaction  and  diffusion) 

Here  flows  1  and  forces  X  in  the  thermodynamics  sense  of  lrrevcrs 
io  <j  processes  will  oe  q^,  7t  ,  r,  J^,  and  corresponding ly 

1  oT  0vt  _  d  r 

T  Os/  ~H7k‘  ~  Eo  —  ^  07x  ~T 

satisfying  the  relation 

o  =  |  Z  IX  (1.11) 

By  assuming  a  linear  relation  of  flows  to  forces,  for  the  iso¬ 
tropic  medium  wc  can  write 


/  Ov. 


,  Or, 

r  —  —  h'n-.  -|-  L--  , 


i  OT 

0  P.; 

fo7~ 

DT 

07, -f 

o 

Oiy 

-  j6lt 

u7, 

}  -j-  C6. 

1  OT 

n  = 

— 

CTo7{ 

Ov, 


(1.12) 


Os,  T 

'  {* 


Here  A,  B,  C,  D,  K',  L,  M  are  kinetic  coefficients  which,  gener¬ 
ally  speaking,  are  functions  of  the  state  of  the  medium,  q  and  £  are 
the  shear  and  second  coefficients  of  viscosity. 

The  requirements  of  Curie's  theorem  were  taken  into  account. 

For  example,  the  tensor  force  dv-j/clx'k  and  scalar  force  e-  cannot  enter 
linearly  into  the  expressions  for  q  and  J,  since  in  both  cases  a 
vector  kinetic  coefficient  would  be  required,  which  is  precluded  for 
the  case  of  an  isotropic  medium.  For. the  same  reasons  the  vector 


force:!  do  not  enter  lntu  .lie  «xprt.-ss  1  on.:  for  v  and  r. 

By  virtue  of  the  principle  of  symmetry  of  kinetic  coefficients 

11  C,  L  -M  (l.li) 

A.-:  .la  apparent,  .In  the  <-xpre;;;;ion  for  there  i;;  a  term  that 

is  proj  ort  lonal  to  e-  and  not  connected  directly  wltli  '.lie  deformation 

of  the  velocity  field,  Thun,  tin-  ..otal  tenaor  of  ntrenne.;  In  the 

relaxing  medium  la  compoaed  of  throe  parta:  l)  "Lhcrmodynami  <• 

prcaauro"  (with  the  ml  nun  sign)  q  -  6  -•  -p-  e  6  j  l?)  vlacoua 

atreanen,  and  y)  the  term  M6^e;.  Jn  accordance  with  (l.ly)  the 

term  Ldv^/ox^  appeara  in  the  expression  for  scalar  flow  _r. 

We  note  that  there  la  no  due  clarity  concerning  the  queation 

of  theae  terms  in  the  literature.  For  example,  De  Groot'n  atudy  [bj 

asserts  „hat  the  velocity  oi‘  the  chemical  reaction  does  not  depend 

on  the  scalar  combination  of  viscosity  forces  dv^/ox^,  since  hydr  - 

static  pressure  does  not  depend  on  the  occurrence  of  a  chemical 

reaction.  It  is  apparent  from  Eqs.  (1.12)  and  (l.ljj)  derived  above 

that  within  the  scope  of  general  phenomenological  considerations 

these  assertions  will  not  be  valid.  In  phenomenological  theory, 

the  question  of  the  magnitude  of  the  aforementioned  effects  remains 

open.  They  can  be  estimated  in  the  molecular-kinetic  theory. 

The  first  two  eq^a iJs.cis  of  (1.2),  Eqs.  (1 .7)  and  (1.8)  together 

with  (1.12),  (l.ly)  form  a  closed  system  of  relaxation  hydrodynamic 

equations.  Here  we  take  p,  s,  and  ^  as  Independent  variables  defining 

the  thermodynamic  state  of  the  medium.  The  remaining  thermodynamic 
\ 

variables  are  determined  from  the  characteristic  function  e(p,  3,  £), 
for  example  T  =  (Be/Bs)p,  . 

When  it  is  possible  to  disregard  viscosity,  thermal  conductivity, 
and  flow  j,  the  system  of  equations  takes  the  form 
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Ji  +  piliv  V  =  0, 
ds  e.i 


pj?  ~  fra(l  P  -  0 
rft* 

PJ7  =  “  A  1 


In  these  conditions  L  =  M  =  0,  which  foil  own,  a::  in  •.•vliicn-  from 


the  condition 


^  A,  A',  A;  >  o 


(1.1b) 


From  Eq.  ( 1 . 1J* )  a  relation  1s  established  between  coolTicient 
K1  and  relaxation  time  t  of  .snail  deviations  from  local  thermodynamic 
equilibrium.  If  in  the  expansion  of  the  right-hand  side  of  the  last 
of  Eqs.  (1.1*0  in  powers  of  (£  -  £o),  we  limit  ourselves  to  linear 
terms  and  if  we  take  into  account  that  the  changes  in  entropy  based 
on  (l.l4)  will  be  of  the  magnitude  of  the  second  order  of  smallness, 
we  obtain 


T  “  AV-«  • 


A'  -  ? 


(1.1b) 


2.  Certain  properties  of  the  equations  of  relaxation  hydrodynamics 
Calculating  the  curl  of  the  right-  and  left-hand  sides  of  the  second 
of  Eqs.  (1.2),  we  obtain 


~  -rourl  (K  X  v)  =  ~  TJP  X  VP  -r  yourl  1*  -  p  VP  x  p 
R  =  curl  V,  P  =  V’h 


(2.1) 


If  we  disregard  viscosity  forces  and  consider  that  hydrostatic 
pressure  depends  only  on  density,  then  (2.1)  transforms  into  the 
familiar  equation  for  an  ideal  fluid 


-fourl  (R  Xv)  -=  0 


(2.2) 


from  which  the  maintenance  of  circulation  folLows.  In  the  presence 


oi'  the  relaxation  pr ocenaoa ,  p  p(p,  a,  .-),  and  lluus  circulation 
la  not  maintained  even  in  tlie  abaence  oi‘  vLacoua  atreaaea  (P  =  0). 

Let  ua  examine  amall  deviation.;  from  the  atate  of  equilibrium. 

By  retaining  in  (2.1)  only  Mu  term.;  of  the  flrat  order  of  .ana  1 .1  no.;.; , 
we  obtain 

■ir  i'  . 

,.  ' .  .  It  IC  curl  ‘ . 

where  v'  la  the  disturbance.;  a'  velocity,  p()  and  q  are  the  undisturbed 
valuea  of  densi.y  and  ahear  vlaenalty. 

Till.;  equation  doea  not  contain  aourcea  and  aincc  the  dlaturbancea 
are  amall  everywhere,  the  flow  in  thla  approximation  la  irrotational . 
By  keeping  in  (2.1)  .he  terma  of  tin.-  aeeond  order  of  amallneaa,  we 
obtain 

<-it*  i  i 

~r  “  r,„5',rl  p  ij  v.r/  v'i-  vf'*  >•  p1  {?-.:) 

In  .he  laat  equation,  the  va.uea  of  the  flrat  order  of  ama  linear 
are  denoted  by  a  aingle  prime  ai  ;n,  while  thoae  of  the  aeconu  order 
of  amallneaa  have  double  prime  aign .  Since  entropy  change a  are  of 
the  aeeond  order  of  amallneaa,  we  can  asaume  that 


VV  =  P>V p‘  -i- 


(2,4) 


Due  to  the  irrotational  character  of  flow  in  the  fir at  approxi¬ 
mation,  we  have 

£v'  =  v(V-v')  (2.3) 

Therefore,  by  net ting  M  ~  0  in  (1.12),  we  obtain 


P'  --  n-dv'  +  ~  nvi.v  - v') 

="  (;(  q-r  C)  V(vv')  - 


-i-  bV(V'v')  •= 


(2.6) 
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When  calculating  P"  we  shall.  consider  fur  simplicity  -he  coeffi¬ 
cients  of  viscosity  lo  be  constant  and  equal  to  their  values  in  an 
unperturbed  medium, 

Then 

1**-  sjAV-i-  l)v(vv")  (2.7) 


Thus 

curl  l'*  —  l)AIC 


Using  (2.3)  and  (2./),  we  write  (2.2)  as: 


Ut 


"si  • 
fV 


■VP'  >:  v„r 


(2.d) 


(2.9) 


It  is  known  that  In  an  ideal  fluid,  disturbances  in  the  velocity 
curl  are  not  displaced  relative  to  the  fluid  (absence  f  transverse 
waves ) . 

This  is  readily  seen  if  we  rewrite  (2.2)  as 

°i  3-  (v-v) 1!  ~  (h-v)v  H- 1!  div  v  =--  0 

For  ohe  case  of  a  purely  rotational  field  div  v  =  0  and  dR'/dt  = 
=  -(R  )v.  It  then  follows  that  disturbances  in  the  velocity  curl 
R  are  "frozen"  into  the  fluid.  The  presence  of  shear  viscosity 
leads  to  a  displacement  of  disturbances  in  the  curl  relative  to  the 
fluid;  this  process,  however,  is  analogous  to  diffusion  and  not  to 
the  propagation  of  waves.  The  right-hand  side  of  (2.9)  defines  the 
speed  with  which  the  eddy  emerges  at  a  given  point  due  both  to 
relaxation  processes  and  viscosity.  In  this  manner,  both  the  relaxa¬ 
tion.  processes  and  the  second  (dllatational)  result  in  an  eddy,  however 
they  cannot  cause  diffusion  of  the  eddy  through  the  fluid. 
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Let  us  consider  the  equations  which  determine  the  propagation 
of  longitudinal  waves  in  a  relaxing  medium.  As  we  have  already  noted, 
in  the  absence  ol'  viscosity  and  thermal  conductivity  the  equations 
01'  relaxation  hydrodynamics  can  be  written  in  the  form  of  (l.l4), 

On  examining  the  small  deviations  from  the  state  of  equilibrium 
and  limiting  ourselves  to  terms  of  the  first  order  of  smailness, 
we  obtain 


Ae;  =  A'k^  1'  -r  At^p', 


ov 


iff. 

0\dft 


(2.10) 


Taking  into  consideration  (2.4)  and  (2.10),  we  write  Eq.  (l.l4) 


as : 


,1,  +  Po  ,77  -  Po  77.  T  a-  7-  Pi.  0  ----  0 


at  1  "•<)/* 


01 


op 


Oi 

i  47, 


l 

iii 


■4(-fr-o  ( 


(2.11) 


-  -  r»:  p. 


Here  ck  is  the  propagation  velocity  of  the  sound  oscillations 
having  a  frequency  u>  »  1/t,  The  primes  on  the  variables  of  the 
first  order  of  smallness  (u,  £,  p)  are  omitted. 

The  formal  solution  of  the  last  equation  of  (2.11)  is  represented 


as 


(2.12) 


Here  A-1  is  an  operator  inverse  to  the  operator  A.  Using  (2.12) 
we  eliminate  £  from  the  second  equation  of  (2.1l)  and  then  taking 
into  consideration  that  operators  A  and  A”1  are  commutative  with 
b/bx^,  we  obtain 


i)vi 

Po  '-Ji  1“  Ccc 


,£p 

dx: 


—  Pi 


—  A*1  — 


It 


()£■ 


=  0 


(2.13) 
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Differentiating  the  first  of  Eqs.  (2.11)  with  respect  to  x 
and  Eq.  (2.1>)  with  respect  to  t_,  we  can  eliminate  p  from  (<•:.  1}). 

C;„  i’ll 


.  tu_Zi  .  -  jJ^L\  ,  .  «u  oh'»  . 

‘  \  ot‘~  ’’  c“  Oi^zJ  +  f  l  eu  dztd,k  ~  0 


(2 .  lJt ) 


Aa  noted  earlier,  in  the  approximation  under  conaideratl on 
the  motion  can  be  considered  irrotatlonal,  and  therefore 


0‘V ,. 


■  = 


(2.15) 


Tl’.e  equilibrium  velocity  of  sound  c0  la  related  to  c  by  the 
relationship  [Jl]: 

<v  -  cj  -  /)=  ( 2  •  1° ) 

u 

Using  (2.15)  and  (2.16),  we  rewrite  ( 2 . 1J1 )  as: 

Equation  (2.1/)  describes  the  propagation  of  longitudinal 
disturbances  in  the  relaxing  medium.  As  t ->  0  and  t ->  co,  Eq.  (2.17) 
transforms  into  the  ordinary  wave  equations  of  acoustics,  in  which 
the  role  of  sound  velocity  is  played  by  c0  and  c^  respectively. 

This  equation  was  obtained  and  analyzed  earlier  [5,  7). 

It  is  of  interest  to  consider  the  characteristics  and  character¬ 
istic  form  of  Eqs.  (l.l4).  If  the  disturbances  are  not  assumed  to 
be  small,  then 


d± 

d  x. 


2  i  3s  ,  0 E 

'  -c«i  V  +  Pn - p  Pc  ~ 

;  “  dxi  1  1  ‘  0xi  '  '  <■  dxi 


dp 


„  dh 


(is  P  dsdp 


(2.18) 


By  using  customary  methods  [8],  we  find  that  the  equations  of 
the  characteristics  for  a  one-dimensional  flow  have  the  following  form 


(lx 

Ti 


—  V  ±Cee, 


dx 

di~  V 


(2.19) 
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Thus,  through  over j  point:  on  the  surface  x,  t  paaa  three 
characteristics,  one  oi’  which  coincldoa  with  the  trajectory  of 
the  particle,  and  the  other  two  coinciding  with  the  "trajectories" 
of  high-frequency  diaturbaneea  (m  »  l/x)  propagating  in  the 
moving  medium.  Introducing  the  operator 


u ,  .1  i>  _  J  i) 

f)l  01  'irl'v±  f“)  ol  ~  .77  r“  o't 


(2.20) 


we  write  -he  relaxation  hydrodynamic  equation  (1.1*1)  in  the  ehurac- 
lerinUc  form 


„  ^ ,  r*  &  .v  .  s  i)  £  it  -  p 

c„-  -bT  ±pea  ln  /-.  ~ur  -r  Pk  --  ■---  x  A‘et  ( -  />< 


A'  . 


'■=  —  A'k; 


(2.21) 


y.  Tranafer  to  equilibrium  hydrodynamics.  Equations  (P.l'j) 
do  net  contain  x.  Therefore,  the  characteristics  of  relaxation 
hydrodynamics  do  not  change  as  x  fends  to  zero  and  do  not  convert 
to  -he  characteristics  of  equilibrium  hydrodynamics  definable  by 
the  velocity  cQ.  Thus,  the  change  to  equilibrium  hydrodynamics 
describing  movement  of  the  medium  for  small  values  of  x,  requires 
special  consideration. 

We  shall  introduce  the  parameter  p  =  cx/L,  where  c_  is  the  veloc¬ 
ity  of  sound  (c0  or  c  ),  L  is  the  characteristic  macroscopic  length 
(the  distance  over  which  the  hydrodynamic  values  markedly  change), 
and,  assuming  p  «  1,  we  shall  consider  the  solution  of  system  (l.l4) 
of  the  form  y  (p.,  r,  t;  p),  where  is  some  kind  of  magnitude  of 
p,  v  etc.  Thus  we  are  considering  flows  in  which  spatial  gradients 
of  the  n-th  order  are  values  of  the  n-th  order  of  smallness. 
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Expanding  -j_  In  a  power  series  with  reapect  to  u 

y  <j  (|ir.  /;  |i)  >/  (|ir.  /)  -i-  p >/  (pr,  /)  -i- . . .  (2.1) 

we  derive  from  (l.l't)  that 

vi  0i  w -  -Af;  (2,2) 

The  Jaac  equation  describes  the  relaxation  proceaa  wiih  charac¬ 
teristic  time  t.  Changes  in  density  and  velocity  field,;,  etc.,  occur 
in  the  second  approximation  which  is  divergent  with  time.  Therefore, 
the  expansion  (y.l)  is  advantageously  used  for  describing  movement 
only  into  a  i\ow  of  sufficiently  small  time  intervals.  The  charac¬ 
teristics  of  the  thus  obtained  equations  are  determined  by  velocity 
c  .  Consequently,  by  means  of  expansion  (2.1)  we  cannot  arrive  at 
the  equations  of  equilibrium  hydrodynamics.  The  difficulties  arising 
here  can  be  due  to  the  small-valued  parameter  (t)  being  .he  coeffi¬ 
cient  for  higher  derivatives  (see,  for  example,  Eq.  (2.1/)). 

To  change  to  equations  of  equilibrium  hydrodynamics,  we  shall 
consider  a  flow  from  the  moment  that  the  relaxation  process  in  the 
zero  approximation  described  by  Eq.  (2.2)  is  virtually  complete 
and  at  each  later  moment,  the  deviations  of  parameter  £  from  equilib¬ 
rium  values  are  determined  by  the  fields  p,  s,  v.  Under  these 

conditions  we  can  write  the  following  expansion  for 

i  =  £»  (p,  *)  -r  nS"1  (|ir;  s,  p,  v)  -f  (pr;  p,  s,  v)  -f  .  .  •  (2.2) 

where  are  functionals  of  the  fields  p(r),  s(r),  v(r).  Thus, 

the  coefficients  of  expansion  (2.2)  depend  on  the  time  Implicitly 
in  terms  of  the  fields  p,  s,  v. 

Under  such  an  assumption  of  the  character  of  the  change  of  £ 

with  time,  it  is  evident  that  the  following  expansion  for  ds/dt  etc., 
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la  alao  valid. 


5,‘  -  M-v‘“  (r:  h 

v)  -r  |»3.V13'  (r;  p.  v)  -  .  .  . 

( 3 . 4 ) 

•,)  -  |.V-»  (r;  P.  .v 

.  v)  pJ\  'i,(r;  p,  s,  v) 

(3.3) 

"1  -  |*A‘‘“(r;  p.  *, 

,  v)  (r;  p,  s,  v)  i  ... 

Ub) 

where 

th<‘  eoef f ielents  In  powers 

of  |i  are  also  functional.: 

of  field; 

p(r). 

s(r),  y(r).  The  explicit 

form  of  these  coefficients 

can  be 

obtained  by  equating  (3.4)  -  (3. 

b)  with  the  expressions  fur 

Op/  v)t 

etc., 

obtainable  from  (l.l4).  For  example,  considering  the. 

t  for 

grad  p  and  dlv  y  there  are  valuea  of  the  order  of  a,  we  find  from 


(3.4)  and  (l.l4) 

|i.Vu)  =  -  (v-grad  p)  -  p  div  v,  .Y1SI  =  A’131  =  .  ,  .  >_  u  (3.7 ) 

Taking  into  conaideration  the  expan3ion  for  p 

/»  =  A».(P.*)-r  M  i"1 -r  ...  (3.3) 

where  p0  la  equilibrium  pressure,  from  (3.l4)  and  (l.l4)  we  find 

Tv‘"  =  -  (v,  V)  v  -  ~  grad  pt,  |xV12’  =  -  j  grad  (3.9) 

From  (3.4)  and  (l.l4)  we  find 

pfin  =  —  (v,  grad  s)  (3.10) 


Leaving  terms  of  the  first  order  with  respect  to  p.  in  expansions 
(3.4)-(3.6)  and  taking  into  account  (3.7)>  (3*10),  we  obtain  the 
equations  of  motion  of  an  ideal  fluid  with  an  equilibrium  value 
of  pressure. 

Note  that  the  change  to  equilibrium  hydrodynamic  equations  by 
expansions  of  type  (3.3)>  (3.4)-(3.6)  are  completely  analogous  to 
the  methods  worked  out  in  nonlinear  mechanics  and  statistical  physics 
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in  connection  with  the  appearance  of  secular  terms  in  expansions  of 
type  (3.1). 


To  obtain  equations  of  motion  In  the  next  approximation,  we 

(0 


need  to  find  '  entering  into  (3.9). 

From  (3.3)-(3.6)  we  find,  correct  to 


a ,,  ih  .\  in  -t-  u^s"' 

ut  “  o;i  '  os 


Substituting  (3.7)  and  (3.10)  into  (3.11),  we  obtain 


v-gr.’ul  to  - orP  <hv  v 


dt 


(3.U) 


(3.12) 


On  the  other  hand,  with  the  same  approximation 


-A’£t=  - 117  5'" 


(3.13) 


Substituting  (3.12)  and  (3-13)  into  the  last  of  Eqs.  (2,21),  we 


determine  % 


(0 


'^‘»  =  £viiw 


By  substituting  this  value  of  into  (3-9),  we  obtain  the 

equation  of  motion  as 


</» 


l>i,77  =  “  Srnd  Po  +  gradftdivv) 


(3.14) 


where 


r  _  _  (0p\  ato 

T  dF  P  whe"  £  =  Tp  (c„=  -  , 


(3.15) 


plays  the  role  of  the  second  coefficient  of  viscosity. 

This  conclusion  assumes  p  «  1.  The  conclusion,  however,  is 
not  restricted  by  the  conditions  of  smallness  of  the  change  in  the 
hydrodynamic  values  and  is  thus  also  valid  outside  the  framework  of 
accoustic  approximation  (compare  the  studies  of  V.  Finkel'nburg  [l] 
and  V.  Granovskiy  [4]), 


The*  propagation  of  small  disturbances  in  relaxing  media  have 
a  number  of  Important  features.  Due  to  dispersion,  a  relaxing  medium 
does  not  have  a  definite  velocity  of  sound,  although  its  properties 
can  be  characterized  by  two  magnitudes  having  a  dimensionality  of 
velocity:  e  and  c0.  The  complete  set  of  these  two  values,  as  will 
oe  shown  below,  l'uiiy  define  the  law  of  propagation  of  .anal!  disturb¬ 
ances  and  thus,  to  a  certain  degree,  replaces  the  velocity  of  sound 
in  relaxation  hydrodynamics. 

In  equilibrium  hydrodynamics ,  small  disturbances  are  propagated 
along  the  characteristics.  In  relaxing  media  this  assumption  uues 
not  hold  by  virtue  of  dispersion,  and  thus  the  characteristic  direc¬ 
tions  take  on  a  different  meaning.  This  has  already  been  seen  from 
the  fact  that  the  characteristic  directions  of  the  equations  of 
relaxation  hydrodynamics  for  all  values  of  relaxation  time  x  are 
determined  by  velocity  cK  simultaneously  with  the  propagation  of 
small  disturbances  for  small  values  of  x. 

To  clarify  both  the  laws  of  propagation  of  small  disturbances 
in  relaxing  media  and  the  role  of  the  characteristic  directions, 
we  shall  consider  the  following  problem  [7]. 

4.  Weak  shock  wave  in  a  relaxing  medium.  Let  a  cylindrical 
tube  be  closed  at  one  end  by  a  piston  which,  when  t  <  0,  is  at  rest 
at  point  x  =  0,  while  at  time  t  =  0  it  is  set  into  motion  and  subse¬ 
quently  proceeds  with  constant  velocity  uQ  in  a  positive  direction 
of  the  x-axls.  The  velocity  of  the  piston  is  small  in  comparison 
with  the  velocity  of  sound  (uQ  <  cQ  <  cM).  Let  us  consider  the 
propagation  of  disturbances  arising  from  piston  motion. 

The  equation  describing  propagation  of  a  small  disturbance  of 
velocity  v  in  one  dimension  has  the  following  form  (see  (2.17)) 
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f  (£.  (Ola  _  4.  c08  Pi  =*o 

T  01  \di>  m  Z?  1  ~  011  0  Otl 


(4.1) 


The  boundary  and  initial  conditions  in  this  problem  can  be 
stated  in  the  following  manner: 


!>(/,  r)-0»!u>n  I  <  o,  y(0,  0  =  /(0  =  {  I!„w"n  (; 


(4.2) 


A  solution  is  found  only  for  the  region  x  >  0.  The  discontinu¬ 
ous  function  f(t)  from  (4.2)  can  be  presented  as  a  contour  integral 


('1 .3) 


Then  the  solution  of  the  problem  is  also  represented  by  a 


contour  integral 


>  (X.  /)  =  -  ^ ^ c/o>,  k  =  0, 

^  <V  —  ic"  (.IT 


(4.4) 


Integration  in  (4. 5)  is  carried  out  along  the  actual  axis  with 
a  by-pass  of  the  origin  of  the  coordinates  with  respect  to  the 
upper  half-plane. 


In  the  lower  half-plane,  function  k(u>)  contains  two  branching 


points 


"l  “  “  ‘  “2  T  *  CD,  =  —  d- 


The  section  for  transferring  to  the  other  sheet  of  the  RIemann 
surface  extends  along  an  Imaginary  axis  between  these  points.  Since 
the  path  of  integration  does  not  pass  through  the  section,  we  can 
select  for  the  entire  path  one  definite  branch  of  the  function  K(co), 
namely  that  for  which  k  =  -kd/c  when  to  ->  0. 
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I 


Lit 


l-’lg.  1. 


It  la  readily  aeon  that  (4.4)  aatlafiea 
both  differential  equation  (4.1)  and  boundary 
condition  (*1.2)  if  we  repreaont  f(t)  in  the 
form  (4,p).  To  verify  fulfillment  of  initial 
condition  (4.2),  we  note  that  in  the  upper- 
part  of  the  half-plane  the  integrand  doou  not 
contain  the  characteristic  features,  then  the 
contour  of  integration  can  be  carried  out  to 
infinity  along  the  positive  direction  if  the 
imaginary  axis  (Fig.  l). 

In  thlo  case 


^  2  Thus,  if  t  <  x/c^,  the  integrand,  when 

integration  is  carried  out  to  infinity, 
contains  an  exponentially  damping  coefficient,  and,  consequently, 
v(x,  t)  =  0  when  t  <  x/c  .  Hence,  in  particular,  the  initial  condi¬ 
tion  (*1.2)  follows  for  x  >  0 

Let  us  consider  the  behavior  of  Jie  solution  along  the  straight 

line  x  -  c  t  =  0,  which  coincides  with  one  of  the  characteristics  of 
00 

Eqs .  (4.1).  As  shown  above;  on  one  side  of  this  straight  line,  namely 
where  t  <  x/c  ,  the  disturbance  is  absent. 

Let  us  consider  the  nature  of  the  disturbance  when  t  -  x/c  = 

=  t1  >  0.  We  deform  the  contour  of  integration  into  a  semi-circle 
of  large  radius  located  in  the  upper  half-plane  with  segments  of  the 
real  axis  adjacent  to  it.  Such  an  Integral,  taken  from  a  mirror 
image  of  this  curve  in  the  lower  half-plane,  equals  zero,  since  on 
Integrating  to  -»  along  the  imaginary  axis  the  Integrand,  when  t1  >  0, 
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decreases  exponentially.  Combining  the  Integrals  with  respect  to 
these  curves,  we  represent  v(x,  t)  as  an  integral  over  a  circle  of 
large  radius  with  its  center  at  the  origin.  In  this  circle  the 
function  k(oi)  can  be  expanded  into  a  series  in  powers  of  (u,t)-1 


k  =  —  { 1 


1  s'_ 

2  iu)t 


1  .1  .\  :.=■ 

2  s' 


)«&+■•■ 


Thus,  if  we  limit  ourselves  to  the  first  three  terms  of  expan¬ 
sion  (4.5),  then 


where 


‘-"-t  h.?) 

By  substituting  variables  w  t'/X  =  X'  ,  we  find 

'  =  (4.8) 

If  we  use  the  familiar  expression  for  Bessel's  function  with 
a  whole  number  index 


Jn  ^  =  liHi  $  CA'P  {y  2  (X  -  t)} 


(M) 


where  the  integration  is  carried  out  along  any  closed  contour  around 
X  =  0,  and  we  rewrite  (4,8)  as 


v(x,  t)  =  a0exp  ^  x'jJt  (—  2  i 


(4.10) 
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To  clarify  the  criteria  of  the  applicability  of  the  obtained 


expression,  we  retain  terms  of  higher  order  in  expansion  ( h .  . 

Aa  a  reault  we  obtain 

i'  (i.  0  =  u„oxp^  -  -r j  |^o  (-  I  *0  i  l  l  Ji  (~  I  >-'  )J  (!|  •  11  ) 


wher 


,1 


The  .n-,:ond  term  in  (4.11)  can  be  negl ected,  if  t 1  a/ A  ~  t. 
Fnerefore  (4.iu)  really  represent:!  the  solution  in  a  narrow  flelu 
(leaa  than  x  in  width)  along  the  characteristic  in  question,  .iineo 
J0(21  yXt 1 )  «=  1  when  t1  «  t,  the  results  thus  obtained  can  be 
written  aa: 


I  n  whon  i'  -  t  —  i  vm<  0 
1  U'.e"'* ,;hon  r>o,  c/(<i 


i 

i 


(Jl  .12) 


Thus,  along  the  characteristic  in  question,  there  arises  a 
discontinuity,  whose  Intensity  decreases  exponentially  with  increas¬ 
ing  distance  from  the  piston.  At  the  distance  of  x  »  c^x  the 
discontinuity  on  the  characteristic  virtually  vanishes  along  with 
the  disturbance  propagating  along  the  characteristic. 

To  obtain  a  complete  representation  of  the  flow  when  x  »  c^x, 
we  must  investigate  the  region  far  from  the  characteristic  (i.e., 
when  t'  »  r).  We  shall  shift  the  path  of  integration  toward  nega¬ 
tive  values  of  the  Imaginary  axis  (Fig.  2).  When  t '  >  0,  the  inte¬ 
gral  with  respect  to  the  path  displaced  toward  -<»  vanishes,  and  the 
solution  is  determined  only  by  the  integrals  around  special  points. 
The  integrand  contains  the  coefficient  which  for  a>  =  aq, 
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('U.1) 


Thu:;  when  x  ;  •  c  t,  the  dlaeont  Inul ty  with  respect  to  ;.}.<■ 
characteristic  gives  way  to  a  continuous  disturbance  whl  *h  arise.; 
along  '.he  straight  line  x  -  el  =  0,  and  whicli  due.;  not  .■  1;.  lm- 
w  1  tii  the  characteristic ,  Tnc  width  of  Lne  region  In  will  m  ul.;tnri,- 
anee  appear;;  in  oi'  the  order  of 

s  (''.Id) 

"  i  IV'  J/ 

anu  there  fore  it.;  relative  ulmen.ilon.:  decrease  with  distance  fr -m 
tnc*  pi.;. .on 

7-/%'  ('1.20) 

Thu.:  'ho  .;mal  ]  perturbation  which  arises  at  a  certain  ;  .in 
J .;  propagated  clone  to  it  (x  <  t/ct  )  a  lout!  tlie  Hiarae  ter  1 .; ..  1  • 
passing  through  this  point;  however,  an  distance  increases,  an  ever 
greater  part  oT  the  perturbation  ia  dlaplaced  l'rom  the  character.!.;'..,!  c 
on  which  it  arlaea,  and  finally  when  x  »  c  i/or  the  propagation  of 
the  perturbation  la  determined  by  the  velocity  c  .  Since  there  are 
no  perturballono  having  a  propagation  velocity  greater  than  <:  , 
the  characterlatlc  x  =  c^t  aeparates  the  unperturbed  region.  However, 
due  to  rapid  damping  of  the  high-frequency  componenta  at  large 
distances  compared  to  c^t,  the  boundary  between  the  perturbed  and 
unperturbed  flow  does  not  proceed  along  this  characteristic. 

According  to  the  result  obtained  from  analyzing  stationary 
solutions  in  Waldmann's  work  [9]  (see  also  Pinkel ’burg 1  a,  et  al. 
work  [l]),  a  weak  shock  wave  in  a  relaxing  medium  does  not  have  to 
contain  discontinuities.  Its  width  is  determined  in  S.  P.  D'yakov's 
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work  [ 10 j .  The  results  obtained  in  this  work  illustrate  a  law,  in 
accordance  with  which  a  weak  discontinuity  in  a  relaxing  medium 
disappears.  Note,  however,  a3  will  be  shown  below,  that  this 
conclusion  i3  not  wholly  applicable  to  plane  flows.  Near  the  barrier 
whose  width  is  comparable  with  the  relaxation  length  l,  and  even  in 
the  presence  of  relaxation  processes,  there  may  arise  a  steady 
discontinuity  of  arbitrarily  small  intensity  which,  however,  vanishes 
exponentially  with  distance  from  the  wall. 

The  unrestricted  widening  of  the  region  in  which  the  disturbance 
arises  (30c  (J1.19))  is  associated  with  neglect  of  nonlinear  effects. 
We  know  that  small  nonlinear  effects  accumulate  and  at  sufficiently 
great  distances  lead  to  marked  deviations  from  a  linear  solution. 

For  instance,  during  propagation  of  a  disturbance  arising  from 
harmonic  oscillations,  the  .leading  edge  of  the  wave  front  gradually 
becomes  steeper.  In  the  case  under  consideration  the  same  effects 
ultimately  limit  the  increase  in  width  of  the  disturbance. 

However,  if  the  disturbance  is  sufficiently  weak  and  limited 
in  time,  it  may  disappear  due  to  absorption  before  a  steady-state 
sets  In.  Actually,  if  at  the  moment  t  =  T  the  piston  stopsj  there 
arises  a  disturbance  which  propagates  according  to  laws  described 
above,  but  at  a  distance  Ax  «  c0T  from  it  (taking  into  consideration 
that  c0  ~  c  ).  Due  to  this,  when 

5  =>/ Mo:Ci>T.c »  c0'/- 

both  disturbances  are  superposed  on  and  weaken  each  other.  Therefore, 
a  disturbance  of  duration  T  disappears  at  distances  significantly 
greater  than  c0T2/m-oT. 
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5 .  Flow  past  a  thin  wedge.  We  shall  consider  the  shape  of  a 
plane,  steady  supersonic  l’low  in  a  relaxing  medium  and,  In  parti eu lar , 
the  pronlem  oi'  the  l'low  past  a  thin  wedge  (having,  an  apex  an(;]r  up 
2qi)  oi‘  a  supersonic  flow  at  zero  angle  oi‘  attack  [  11 ,  12],  Lot  he 
x-aels  coincide  with  the  direction  oi’  movement  oh  an  unpcrturoeU 
oW  whose  velocity  we  denote  by  U,  and  let  us  assume  that  1)  )  p  c  ,■ 

>  ■  The  equations  oi‘  relaxation  hydrodynamics  for  the  plane, 

steauy  flow  are  written  as  (see  (l.l^l)): 
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where  v  ,  v  are  the  Velocity  components  a  ion,'  the  x-  and  y-axes. 

>■  j 

in  ..hose  equations,  due  to  thermal  conductivity  and  diffusion,  tne 
anisotropic  portion  of  the  stress  tensor  and  change  in  entropy  can 
no  see  looted.  The  influence  of  these  effects  is  evaluated  below. 

since  the  perturbations  caused  by  the  wedge  are  negligible, 
Eq.  (6.1)  can  be  linearized  by  setting  v  =  U  +  u,  v  =  v,  etc. 

**  t/ 

Thus  we  obtain  for  v 


Here  M  =  u/c  ,  FL  =  u/c.  are  the  Mach  numbers.  Similar  cqua- 
tions  can  be  obtained  for  density  and  velocity  disturbances  along 
the  x-axls . 

Equation  (5.2)  is  of  the  hyperbolic  type  for  Mra  >  1  and  of  the 
elliptical  type  for  M  <  1.  In  the  first  case,  three  characteristics 
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paoa  through  each  point,  two  of  which  have  the  form 

*  -  ±  V^lca*  -  \’J  -I-  const  (  5 ,  '  ) 

while  i he  third  la  parallel  to  the  x-axia.  The  direction  or 
theae  characterlatlca  doea  not  depend  on  t  and  doea  not  coincide  with 
the  direction  of  the  Mach  llnea  of  equilibrium  hydrodynamic;;, 
x  =  -  ^  Mq  -  ly  +  c o n 3 1 .  In  order  that  the  function  exp  (Bx  -  Dy) 

represent  the  partial  solution  (5.2),  it  ia  necessary  that 


n  ‘H-vj-un  -w-i) 

D  dll  ~  i 


(5.'0 


We  will  consider  the  origin  aa  coinciding  with  the  apex  of 
the  wedge  and  we  shall  write  the  boundary  conditions  of  the  problem  aa 


I  0  when  r<0 
V  M  =  {<p U  when  x>0 


(5.5) 


Representing  f(x)  aa  a  contour  integral,  we  will  write  the 


solution  of  the  problem  in  the  following  form: 


v  {'x'  y) =  ^  ^  ~ a  (11)  y>^ dr] 


where 


(5.6) 


y  iii+i 


=  L  w  =  !L 
x  i  •  y  i 


(5.7) 


In  (5.6),  integration  is  carried  out  along  the  real  axis  with 
a  bypass  of  origin  of  the  coordinates  in  the  lower  half-plane. 

The  function  a(q)  in  the  plane  of  the  complex  variable  is 
analogous  according  in  its  own  properties  to  the  function  k(to) 
introduced  in  Section  special  points  of  this  function  are  distrib¬ 
uted  in  the  upper  half-plane  t)i  =  i,  =  l(M^  -  l)/(M^  -  l).  The 
path  of  integration  in  this  case  also  passes  through  the  cross 
section.  When  q  ->  «  and  tj -*  0,  a(q)  takes  on  the  real  values 
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<i  M  =  ±  YMJ  -  i. 


a  (0)  =  t/iW.1  -1 


IT  we  select  that  branch  of  (5.7)  which  a  a  i)  ->  ®  approaches 

t  ^M'  -  1  ,  then  by  transferring  the  contour  of  Into, ration  -» 
alone  he  imaginary  axis,  we  can  ah<»w  that  the  solution  vanish'-,;  wii< -i 

i-f^Ty<0  (L.j) 

o'.  herwlse  vanishing  >.ako:s  place  when 

,•  ■  I  '.Y/.;*-  ty  •  II  (-  • lu) 

Sunaequent ly  in  the  solution  l'or  y  >  0  it  la  nee  canary  j  une 
tiiat  uranc.n  of  the  radical  (5.7)  which  corresponds  to  the  plus  sign 
in.  (5.d) .  The  second  branch  should  yield  perturbations  propagat J  n. : 
upstream  and  which  should  be  discarded  for  physical  oonsideratl  o.ns. 
yhea  'j  <  0,  we  must  choose,  for  the  same  considerations,  a  nranel, 
corresponding  to  the  negative  sign  in  (5-8). 

Equation  (5.9)  is  one  of  the  characteristics  passing  thr  >u.gh 
the  apex  of  the  wedge.  Let  us  consider  the  behavior  of  the  solution 
a.ong  this  characteristic  in  a  layer  of  width  «  1.  Deforming  the 
contour  of  integration  in  the  same  manner  used  in  the  problem  of 
.he  piston  motion,  we  find 


v  (-0  y) 


when 


*  -  K.T/^  -  l//  <  o  (5.H) 


i  VU  CXl)  l-  '4P2  (.'/  I  l))J o  (=')  whnn  *  -  V~UJ  -  1  y  >  0 


Here  J  (z  1 )  is  Bessel's  function 


V-  £^>0. 


'(2,1s 


1  4  K'V-i  r 


-  tv 


The  function  J0(z')  *  1  in  the  vicinity  of  the  characteristic 

(9.9). 
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Thus  near  the  apex  of  the  wedge  along  the  characteristic  passing 
through  it,  a  discontinuity  arlaea,  which,  however,  exponentially 

decreaaea  with  distance,  When  x  »  1,  y  »  l  the  perturbation  along 
the  characteristic  virtually  disappears,  Here  the  main  contribution 
is  determined  by  the  long -wave  components  of  the  perturbation.  If 
in  (b.6)  we  let  the  quantity  a(q)  be  approximately  |4i'  -  1,  then 
after  integrating  we  obtain  a  result  which  follows  from  equilibrium 
hydrodynamics 


I'  (T,  i)  - 


fn 

U -U 


when  1  —  Ma‘  --  1  y  <  0 
when  i'  —  1  ■  1  y  0 


(5.12) 


However  if  in  a(q)  we  take  into  account  the  terms  of  a  higher 
order  with  respect  to  q,  we  find  that  actually  there  will  be  no 
discontinuity  along  the  line  x  -  yM2  -  ly  =  0.  Proceeding  in 
the  same  manner  as  in  the  problem  of  the  piston  motion,  we  find 


Here  3>x  is  the  integral  of  probability  (4.17) 


W  ‘^00*  Q 


(5.15) 


(5.14) 


Equations  (5.15)  and  (5.14)  yield  the  shape  of  the  Mach  line 
produced  by  the  effect  of  relaxation  processes.  The  characteristic 
width  of  the  line  is  of  the  order  of  fv2ly. 

The  flow  under  consideration  can  be  integrated  as  a  weak  shock 
wave  arising  near  the  wedge  in  a  supersonic  flow.  As  noted  earlier 
[9#  lj  in  a  relaxing  medium  there  can  be  no  discontinuities  of  small 
magnitude  propagating  with  a  velocity  less  than  c  to  the  unperturbed 
medium.  In  the  problem  in  question  we  are  concerned  with  the  more 
involved  case  where  the  value  of  perturbation  (cpU)  is  small,  but 


-29- 
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Ju?  velocity  of  the  wedge  la  greater  than  c  .  In  thin  cane,  as 

CO 

evident  from  the  previously  derived  solution  (I3.ll),  in  the  vicinity 
oi  tiic  wall  and  in  the  region  ol'  the  order  o;‘  l  Ut  a  steady 
discontinuity  ol’  arbitrarily  small  Intensity  can  arise.  Within 
oln-  r'ramewurk  ui’  relaxation  hydrodynamics,  we  have  to  consider  tnls 
discontinuity  as  a  geometrical  surl’ace.  Its  width  is  detenu.! nod 
ay  tne  Inie  oi  the  Tree-path  length. 

Another  feature  oi  these  discontinuities  which  arise  near  the 
narrier  is  that  the  discontinuities  are  directed  (lor  small  Intensi¬ 
ties)  along  the  characteristics  of  the  relaxation  hydrodynamic 
equations  and  thus  do  not  coincide  with  ordinary  Mach  lines.  Yhey 
vanish  when  the  valoclty  of  the  wedge  becomes  less  than  c  ,  but  not 
less  than  u(j.  This  is  because  Eq.  (5-2)  becomes  elliptical  when 

U  -1  :  . 

00 

To  elicit  the  possibility  ol'  experimental  detection  ol  those 
discontinuities,  we  must  evaluate  the  said  assumptions.  Basically 
they  reduce  to  neglect  of  the  effects  of  viscosity  and  thermal 
conductivity.  When  the  solution  contains  frequencies  close  to 
relaxation,  we  can  neglect  the  dissipation  effects  associated  with 
viscosity  and  thermal  conductivity.  Actually,  the  effect  exerted  by 
the  latter  Is  of  the  order  of  Zq/l  =  tqc/l  (lQ  is  free-path  length, 
t  ^  Is  free-path  length  time,  £  is  average  thermal  velocity  of  the 
molecules,  L  is  the  characteristic  dimension  of  the  region  In  question), 
on  the  other  hand,  oh o  influence  of  the  relaxation  effects  is  of 
the  order  of  l/h  =  Ut/l  »  i  /L,  since  t  »  tq. 

With  distance  from  the  wall  the  Intensity  of  the  discontinuities 
under  consideration  decrease  exponentially,  and  thus  they  can  be 
detected  experimentally  only  provided  the  "relaxation  length". 
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I  =  ( cm/ c q )  U  t,  la  significant .  In  particular  it  la  necessary 
that  l  be  significantly  greater  than  the  thickness  or  the  boundary 
layer.  We  readily  aee  that  under  normal  preaaure,  the  thlckneaa 
ol'  the  boundary  layer  generally  exceeda  the  dlmenaiona  oi’  the 
relaxation  boundary.  We  know,  however,  that  the  relationship  oi' 
uhe  boundary  layer  A  to  the  characteristic  dlmenaiona  ol'  the  stream¬ 
lined  body  L  equal  a 


where  R  ia  the  Reynolda  number.  Since  q  ia  approximately  l/'j>  clQ, 
A  «  lQL/iM.  Since  the  relatlonahip  l /l  =  N  »  1  ia  virtually 
independent  of  denalty,  the  relationship 


\ 

/ 


„i/  ri  i  f  'L 
>  ;u/nv  V  i„ 


decreases  aa  the  free-path  length  increases,  i.e.,  the  dimensions 
of  the  relaxation  zone  increase  on  rarefaction  more  rapidly  than 
the  thickness  of  the  boundary  layer.  Here  it  i3  absolutely  necessary 
that  the  free-path  length  remain  small  in  comparison  with  L.  Because 
l  usually  exceeds  lQ  by  several  orders,  l Jl  =  l/N  is  a  small  value 
and  thus  the  conditions  are  possible  when 

a  ==  _  L- 

1  A'  V  /0^ 


despite  the  fact  that  L /l  »  1. 

In  the  case  of  strong  shock  waves,  the  behavior  of  the  shock 
wave  in  a  relaxation  zone  near  the  wall  should  be  specially  investi¬ 
gated.  From  similarity  considerations  we  can  expect  that  in  the 
vicinity  of  the  wedge  at  distance  x  «  l,  the  flow  should  be  analogous 
to  a  gas  flow  having  a  frozen  relaxation  process.  Thus  there  should 
arise  a  discontinuity  whose  direction  is  determined  by  the  Mach  number 
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Mm  =  U/c^.  At  the  name  time,  at  distances  x  »  l  I'rom  the  apex 
there  should  occur  an  ordinary  a  hock  wave  (expunde  :  e.e  1 :  the 

effect  oi  relaxation  processes),  whoso  angle  of  inclination  la 
determined  by  the  velocity  <jT  sound  c  .  Since  the  lending  front 
oi’  -hlc.  „ype  wave  has  a  discontinuity  whose  width  i.;  r  tiie  order 
or  the  rree-path  length  [9,  l],  we  can  also  expect  that  In  the 
case  o r  a  strong  shock  wave  the  leading  front  is  curved  as  1; 
approaches  the  apex  of  the  wedge,  gradually  assuming  a  direction 
characteristic  of  a  shock  wave  in  a  medium  lacking  Internal  degrees 
oi“  rreedom.  The  characteristic  dimension  within  whose  limits 
curving  is  noticeable  should  be  of  the  order  of  l. 

In  the  relaxation  zone,  the  disturbances  are  propagated  a.ong 
directions  differing  from  those  which  follow  from  equilibrium 
hydrodynamic  considerations.  At  distances  significantly  exceeding 
the  zone  uf  relaxation,  propagation  is  localized  in  a  relatively 
narrow  region  along  the  characteristics  of  equilibrium  hydrodynamics. 
The  width  of  this  region  increases  proportionally  as  the  square  root 
of  distance,  although  it  doe3  not  begin  to  exert  an  influence  on 
the  nonlinear  effects.  Thus,  the  perturbation  issuing  from  two 
carriers  at  a  distance  L  from  each  other  and  which  merge  into  a 
single  Mach  line  when  y  >  L2/V"l*  In  exactly  the  same  manner  the 
harrier  having  linear  dimensions  L  yields  a  Mach  line  of  finite 
Length,  which  vanishes  at  a  distance  of  the  order  of  L~/v:'l  from 
the  x-axis  due  to  superposition  of  disturbances  issuing  from  extreme 
points.  We  shall  consider  In  greater  detail  the  change  in  the  nature 
of  the  disturbance  as  distance  from  the  site  of  its  origin  increases. 

b .  Propagation  of  small  disturbances  In  the  relaxation  zone. 
Suppose  on  the  line  of  flow,  y  =  0,  there  arises  an  arbitrary  dlsturb- 
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ance  of  some  hydrodynamic  magnitude  (velocity,  density,  temperature) 


/•'  ( i .  0)  -  i|>  U) 


(6.1) 


Ii’  the  di3tir  bance  la  atoady,  it  La  dcacribed  by  the  following 
equation: 


iKJV-ng- 


,  o-y 


n'F 

0\jl 


-  0 


(b.2) 


which  muat  bo  Integrated  wi th  boundary  condition  (u.l).  In  the  caae 
of  equilibrium  hydrodynamica,  the  dlaturbance  la  governed  by  the 
wave  equation 


W- 


1 1  "l! _ 

1 '  Or1  ox/ 


(6.5) 


which  may  be  obtained  from  (6.2),  having  aet  the  relaxation  length 
equal  to  aero.  If  we  discard  disturbances  propagating  upstream, 
the  solution  of  (6.5)  for  y  >  0  is  then  written  as 

F  <*•  y)  =  (*  -  VJQ  -  i y)  (6.4) 

Solution  of  (6.2)  can  be  written  as: 


where 


(6.5) 


(6.6) 


and  function  a(q)  was  taken  from  (5-7).  When  y  »  1,  only  the 
low-frequency  components  of  the  initial  disturbance  f[x)  contribute 
to  (6.5). 

Expanding  a(q)  into  a  series  in  q 

0  Ol)  =  —  IT|  -£  +  ■  •  • 

,  * 

and  limiting  ourselves  to  the  first  two  terms  of  the  expansion,  we 
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I'ind 


f  *  Tn  \\  0XV  I'l  (f  ~  ^  “  1  f  “  f ))  * 

x  exp(-  n’-f-v-)1!1  (OjfVC  = 

t  Cl  I  C  1  V*  t/ I  i|(r  —  |/’A/o*  —  ly. —  C)  i  ir 

■  Ja  JJ  PNi |M1  I  t  TTj  — - L~1 - - — —  d'\ < 


where  i.;  also  determined  In  the  name  way  an  In  (k.1'0 
integrating  with  reapect  to  q,  we  obtain 

f  '  ■  • »>  •  v  k  \  yk,  {-  *£.)♦  <' 1  KST^T,  -  q  </; 


liquation  (o./)  more  exactly  deacribea  the  lawn  of  propagation 
of  small  diaturbancea  for  y  »  1  than  doea  Eq.  (6.4).  Since  in 
oqui linrium  hydrodynamica  the  diaturbance  at  any  point  F(x,  j)  la 
fully  defined  by  the  value  of  the  diaturbance  at  a  definite  point 
P  (x  .  0)  of  the  x-axia,  whoac  coordinate  ia  x  =  x  -  \  M  -  ly. 

If  we  ake  Into  account  the  relaxation  processes,  the  dlaturoance, 
aa  evidouv  from  (o.y),  at  any  point  ia  represented  aa  she  .nun  of 
diaturbancea  tranamitted  by  different  pointa  of  the  x-axia,  however, 
the  weight  with  which  the  contributions  of  the  different  pi inta  are 
introduced,  decrease  exponentially  with  distance  from  the  point 
P  (x(;,  u);  the  width  Ax  of  the  region  of  the  x-axia  influencing 
the  diaturbance  at  point  P  is  of  the  order  of 

Ax  =  ^2lyv2 


It  also  follows  from  (6.7)  that  the  initial  disturbance  f(x) 
is  transmitted  with  a  distortion:  the  integration  results  in  a 
gradual  smoothing  out  of  the  Initial  disturbance  and  disappearance 
of  fine  details  as  distance  from  the  x-axis  Increases. 
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THE  BREAKUP  OE  DROPS  AMD  A  FLUID  STREAK 
BY  AN  A  JR  SHOCK  WAVE 

A.  A.  Euzukov 
(Novsiblrsk) 

Wh!  if  studying  the  disintegration  of  drops  and  a  fluid  stivam  by 
a  sh >ok  wave,  Jt  was  found  [1-2]  that  disintegration  results  from  the 
nr  tracted  effect  of  the  dynamic  pressure  of  gas  flow  behind  a  shock¬ 
wave  fruit,  and  that  the  relative  velocity  of  gas  flow  and  drop  size 
considerably  affects  the  nature  of  disintegration.  Here  the  basic 
defining  parameter  is  the  relationship  between  the  pressure  forces  of 
gas  fl-ws  on  the  drop  and  the  surface  tension  forces,  equal  to  pu2d/o. 
Here  u  is  the  relative  velocity  of  gas  flow,  p  is  gas  density,  d  is 
the  parameter  of  the  liquid  drop,  and  o  is  the  coefficient  of  the 
surface  tension.  From  evidence  in  the  literature  [j-5],  disintegra¬ 
tion  does  not  occur  at  certain  less-than-critical  values  of  this 
parameter,  and  the  liquid  drops  are  stable  to  a  gas  flow  of  subcritical 
velocity.  This  is  because  the  capillary  forces  which  impede  break¬ 
down  of  the  drop's  spherical  stability  predominate  over  the  dynamic 
pressure  straining  forces. 


At  relatively  low  supercritical  flow  velocities,  the  strains 
produced  by  dynamic  pressure  increase  and  the  drop  begins  to  break  up. 
The  mechanism  of  this  breakup  is  as  follows:  the  windward  side  of 
the  drop  Is  pressed  Inward  and  the  drop  becomes  parachute-like.  The 
thus  formed  liquid  film,  convex  In  the  flow  direction,  bursts,  form¬ 
ing  numerous  minute  droplets.  In  certain  cases,  dr  ps  can  be  broken 
up  by  Increasing  oscillations  f  the  drop. 

If  the  velocity  A'  the  gas  fl  >w  greatly  exceeds  critical  velocity, 
the  drop  disintegrates  by  separation  of  the  windward  surface  layer  of 
the  liquid  from  the  drop,  whereas  its  bulk  is  displaced  and  deformed 
negligibly.  It  was  noted  that  this  disintegration  process  Is  preceded 
by  the  emergence  of  minute  capillary  waves  on  the  drop  surface. 

Whereas  the  mechanism  of  drop  breakup  at  critical  and  low  super¬ 
critical  velocities  has  been  sufficiently  well  studied  and  can  be 
theoretically  described  [6-8],  the  circumstances  of  the  occurrence  A' 
"supercritical "  breakup  at  gas  flow  velocities  comparable  to  sound 
velocities  are  insufficiently  clear. 

Described  below  are  the  results  of  an  experimental  study  of  super¬ 
critical  breakup  of  droplets  and  streams  of  water  by  an  air  shock 
wave.  We  also  studied  the  qualitative  picture  of  breakup  and  have 
proposed  certain  empirically  verified  formulas  for  breakup  time.  A 
value  was  found  for  the  critical  relationship  between  the  governing 
parameters  at  which  the  aforementioned  supercritical  mechanism  of 
disintegration  starts. 

1.  Description  of  experimental  setup.  Figure  1  shows  photo¬ 
graphs  of  the  breakup  process  of  water  drops  at  different  instants 
after  passage  of  the  shock  wave  (the  numbers  below  the  picture  Indicate 
in  microseconds  the  time  elapsed  from  the  start  of  shock  wave  flow 
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past  the  drop),  and  FJg.  2  ahows  photographs  of  the  breakup  process 
f  water  streams  by  a  gas  How  behind  the  shock  wave  front:  ( 0 - 1 )  Is 
th<-  time  of  the  preliminary  breakup  stage,  (1-2)  is  the  time  of  the 
actual  breakup.  In  the  photographs,  the  velocity  >f  gas  f!  w  Is 
designated  in  m/sec  (the  number  6  In  the  last  photograph  e  -rresponda 
t  •  a  Vf.  ■  1  t.y  >d’  172  m/sec). 

These  photographs  were  obtained  in  a  shock  tube  (Fig.  3)  f 
•  ns.t.ai.t  :ross  sectl  -n  (110  x  110  mm),  which  ensured  a  greater  than 
2  msec  duration  of  virtually  constant  pressure  and  gas  flow  velocity 
l.n  the  sh'-ck  wave.  This  duration  of  gas  flow  of  constant  parameters 
was  entirely  sufficient  for  completing  the  entire  breakup  process. 

The  gas  fh-w  velocity  behind  the  shock  wave  front,  which  was  varied  In 
different  experiments  from  50  to  200  m/sec,  was  determined  by  the 
pressure  in  the  high-pressure  section,  which  was  separated  l'r  m  the 
working  part  of  the  shock  tube  by  a  cellophane  diaphragm.  T'.  obtain 
the  required  pressure  difference,  the  diaphragm  was  shot  out  of  an 
air  gun.  A  diagram  of  the  experimental  device  is  shown  in  Fig.  3, 
where  1  is  the  working  part  of  the  shock  tube,  2  the  high-pressure 
section,  3  the  water  inlet  nozzle,  4  the  pressure-sensitive  detector, 

5  a  double  pulse  oscillograph,  6  a  delay  line  spark  light  source,  7 
the  aperture  of  the  phototracer,  and  8  the  photorecorder  magazine. 

A  part  of  the  shock  tube  was  made  out  of  plexiglass  in  order  to 
see  the  breakup  process.  The  photographs  were  made  by  the  open 
shutter  method,  using  a  spark  impulse  bias  lighting  triggered  at 
regulatable  time  intervals  from  pressure  detectors.  Moreover,  to 
obtain  a  continuous  time  exposure  of  the  breakup  process  we  used  a 
camera  which  enabled  us  to  evolve  an  image  on  a  film  moving  at  a 
speed  of  100  m/sec. 


Along  with  photographing  the  processes,  shock-wave  pressure 
diagrams  were  oscl 1 lographed  for  a  more  accurate  determination  of 
flow  velocity  and  ga3  density. 

The  droplets  and  streamlets  of  water,  whose  diameter  in  different 
experiments  varied  from  0.8  to  3.5  mm,  were  created  by  U3lng  droppers 
inserted  into  the  shock  tube  through  the  upper  wall. 

2.  Experimental  results.  It  la  apparent  from  the  pictures  in 
Figs.  1  and  2  that  from  the  moment  of  the  passage  of  the  shock  wave 
past  the  drops  or  water  streams,  the  entire  disintegration  process  is 
distinctly  divided  into  two  stages. 

During  the  first  stage,  the  preliminary  stage,  the  drops  or 
water  streams  remain  intact,  however,  disturbances  of  the  surface 
capillary  wave  type  occur  on  their  surface.  These  disturbances 
Increase  with  time  and  are  transformed  into  stable  crests,  whose 
crowns  are  on  the  side  surfaces  of  the  drops  and  streams. 

The  formation  of  these  crests,  from  which  a  gradual  separation 
of  water  begins,  concludes  the  first  stage  —  preparation  for  breakup. 
Depending  on  size,  the  shape  of  deformed  drop  and  the  number  of  crests 
can  differ.  For  example,  while  relatively  large  drops  of  water  with 
diameters  of  3-4  mm  have  two  or  three  series  of  crests  from  each  of 
which  disintegration  of  the  water  occurs,  drops  with  diameters  up  to 
1  mm  are  surrounded  by  only  a  single  crown  of  the  crests. 

The  second  stage  of  breakup  (1-2  in  Fig.  2),  which  lasts  until 
total  disappearance  of  the  drop  or  stream  is  the  stage  of  actual 
disintegration.  During  this  stage,  there  forms  a  jet  of  atomized 
water  consisting  of  almost  imperceptible  droplets  and  vapors.  No 
special  study  of  the  degree  of  atomization  was  made,  but  Individual 
photographs  under  the  microscope  show  that  the  atomization  jet  contains 


liquid  droplets  with  dimensions  of  20-50  p.  During  the  second  stage 
we  observed  nonunl  forml  ty  In  (.he  atomisation  process .  By  the  degree 
■  <;'  opacity  of  the  atomisation  Jet,  It  Is  possible  to  distinguish 
between  the  Initial  ejection  of  water,  evidently  associated  with 
final  formation  I'  crests,  the  Increase  of  flow  rate  t  maximum,  and 
the  final,  raj) Id  diminishing  of  flow  rate  by  the  end  i  atomisation. 

The  times  f  the  first  and  second  stages  depend  essentially  on 
ve .  city  A'  gas  flow  and  diameter  of  water  streams;  when  fl  w  velocity 
drops  t  i  bO-'fO  m/s ec  and  stream  diameter  decreases  to  1  mm,  the  change 
lii  the  mechanism  of  atomization  is  noticeable.  The  breakup  of  the 
streams  begins  to  resemble  disintegration  into  relatively  large  pieces, 
and  not  finely  dispersed  atomization.  An  example  of  this  type  of 
disintegration  is  shown  in  the  fburth  photograph  from  the  top  In  Fig.  2. 

It  must  bo  pointed  out  that  at  the  gas  flow  veJ  cities  being 
smoldered,  both  the  times  of  the  first  and  second  stages  of  atomiza¬ 
tion  are  so  small  that  there  is  insufficient  time  for  the  streams  r 
drops  to  be  substantially  displaced  and  the  internal  streams  do  not 
have  enough  time  to  develop  so  intensely  in  them.  Therefore,  longi¬ 
tudinal  instability  of  the  drop,  which  generally  influences  the  dis¬ 
integration  process,  Is  less  substantial.  In  addition,  this  explains 
the  fact  that  the  flow  rate  proceeds  from  the  windward  side  while  the 
leeward  side  of  the  drops  and  streams  remain  undisturbed  for  a  long 
t  J  me  . 

j>.  Atomization  time  of  drops  and  fluid  streams.  To  calculate 
the  time  of  the  preparation  stage  of  atomization  during  which  surface 
waves  and  crests  form  on  the  surface  of  drops  or  the  stream  we  use 
the  equation  of  steady  motion  of  a  fluid 

<l  U  | 

JT  "t  (» Brail)  u  =  -  j  gr-ad  I>  +  viu 


(3.1) 


Considering  water  as  an  Ideal  liquid,  In  the  equation  wo  can 
disregard  the  term  which  takes  1 n t  -  account  viscous  friction  forces. 
Since  the  surface  crests  are  small  relative  t  the  m;r-:t!l  dlmcna!  i 
uf  the  dr<p  or  stream,  we  can  also  discard  the  second  term  of  the 
equation.  In  this  case,  Eq.  (3.1)  takes  the  form 


,tv  LT-lM* 

dt  '•* 


(3.2) 


Let  us  estimate  the  order  of  the  left- and  right-hand  sides  jf 
Eq.  (3.2).  For1  the  left-hand  side  we  have 


(3.3) 


where  the  diameter  of  the  drop  .>r  stream  d  Is  taken  as  the  character¬ 
istic  dimension,  and  where  the  time  of  the  preparation  stage  of  atom¬ 
ization  ti  is  taken  as  the  characteristic  time.  Since  here  we  are 
considering  a  mechanism  of  atomization  corresponding  to  large  values 
of  the  parameter  pu2d/o,  we  car.  conclude  that  the  chief  force  deform¬ 
ing  the  drop  surface  is  that  of  the  dynamic  pressure  of  the  oncoming 
gas  flow.  By  taking  this  into  account,  Eq.  (3.2)  is  presented  as 

(3.4) 

where  p#  is  the  density  of  water,  p  is  the  density  of  gas,  and  u  is 
the  relative  velocity  of  the  gas  flow.  Combining  (3.2),  (3.3)  and 
(3.4),  we  finally  obtain 

'.-‘.Tj/’j  (3.5) 


An  analogous  dependence  for  the  disintegration  time  nf  a  liquid 
drop  in  a  gas  Is  carried  out  in  Gordon's  work  [6],  where  the  propor¬ 
tionality  factor  is  taken  to  be  2.  As  follows  from  our  measurements, 


thl::  value  for-  the  factor*  la  overestimated .  This  la  understandable 
ali, -i-  tiie  author  In  hla  calculations  assumed  the  deformations  of  tile 
dr  |  de termini ii|.'  the  time  required  for  disintegration  to  be  comparable 
w! tii  the  Initial  diameter,  an  assumption  which  In  the  ease  f  "super¬ 
critical"  atomisation  la  false. 

Ac  ail  wn  empirically,  during  the  time  "f  actual  at,  ,m!.:at  1  m,  the 
.'.quid  I'l.m  pee  i  a  fr  mi  the  windward  aide  of  the  drop  r  stream.  Here 
it  la  appropriate  to  tonalder  the  actual  atomisation  pr  ■cess  as  a 
pr  :esa  f  the  spreading  of  water  over  the  circumfluent  gas  fl jw. 

Tills  spreading  Is  observed  when  the  fluid  is  pierced  by  a  liquid 
"hammerhead."  Therefore,  to  determine  the  actual  atomisation  time 
we  :a:,  use  the  expression  for  the  depth  of  the  opening  being  per¬ 
forated  by  a  cumultative  stream  [9] 

1.  1  \Z  '■  (5.'-) 

Mere  l  Is  length  of  stream,  p*  Is  density  of  the  stream  matter, 
und  p  is  density  of  the  matter  comprising  the  barrier.  Identifying 
the  depth  of  the  opening  with  the  length  of  the  drop  in  the  gas  flow, 
wo  have 

L  ~  lj  u 

(3.7) 


Hence 


*3 


(3.8) 


Here  the  diameter  d_  of  the  drop  or  stream  is  taken  as  the  char¬ 
acteristic  dimension  of  the  drop  or  stream. 

h .  Comparison  with  the  experiment.  For  verifying  the  obtained 
calculated  relationships,  we  carried  out  a  series  of  experiments  to 


determine  the  time  of  the  preliminary  a tage  and  the  actual  atomization 
time  by  measuring  the  diameters  of  the  streams  and  the  gas  flow  veloc¬ 
ities  within  a  wide  range. 

On  photographs  similar  to  those  shown  In  Fig.  2,  the  following 
\ 

measurements  were  carried  out:  the  preliminary  stage  time,  tj.,  was 
determined,  counting  It  from  the  Instant  the  shock  wave  (o)  flows 
past  the  stream  up  to  the  Instant  of  the  appearance  of  the  first 
ejection  of  water  (i);  the  time  of  actual  atomization,  t2,  was  deter¬ 
mined  from  the  appearance  <f  the  first  ejection  of  water  (l)  until 
a  maximum  flow  rate  of  water  or  until  breakup  of  the  stream  Is 
achieved  (2).  In  each  experiment,  we  recorded  the  initial  diameter 
of  the  stream  and  the  velocity  of  the  shock  wave  front,  thus  making 
It  possible  to  calculate  the  proportionality  factors  ki  and  k2  in 
Formulas  (3.5)  and  (3.8). 

The  results  of  measuring  ki  and  k2  in  relation  to  parameter 
pu2d/o  are  shown  in  Figs.  4  and  5.  It  turns  out  that  at  values  of 
this  parameter  of  103  and  greater,  which  correspond  to  flow  velocities 
higher  than  6O-7O  m/sec  and  to  diameters  of  the  water  streams  greater 
than  1  mm,  the  values  of  ki  and  k2  vary  negligibly,  thus  indicating 
the  correct  choice  of  force  factors  used  in  evaluating  the  atomization 
time.  With  decreasing  velocity  of  the  gas  flow  and  with  diameters 
less  than  the  indicated  limits,  the  values  of  ki  and  k  start  to 
increase  sharply,  because  of  the  increasing  effect  of  the  forces  ex¬ 
erted  by  surface  tension  during  deformation  and  breakup  of  water 
streams.  It  is  evident  here  that  the  proposed  mechanism  of  atomiza¬ 
tion  becomes  invalid. 

Experiments  on  spontaneous  measurements  of  atomization  time  of 
sperical  drops  were  not  conducted,  although  the  dimensional  relation- 


nhlpn  (?.5)  and  (3.0)  alao  hold  for  them  with  a 
ic i  and  k(?  on  the  decreasing  aide. 
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THE  STRUCTURE  OF  A  FLAME  FRONT 

A.  S.  Pieshanov 
(Moscow) 


We  found  In  the  first  two  approximations  (from  the  square  of  the 
Mach  number,  M2)  solutions  I'  hydrodynamic  equations  with  respect  to 
chemical  reactions.  We  noted  that  in  the  stationary  case  the  pressure 
in  front  of  the  flame  front,  which  varies  prop  irtJonal ly  as  M2, 
Increases  when  P  >  1  or  decreases  when  P  <  1  (P  is  the  Prandtl  number). 
For  the  nonstationary  case,  qualitative  distributions  of  temperature, 
pressure,  density,  and  velocity  are  given  with  respect  to  the  flow. 

In  the  presence  of  diffusion  processes  and  chemical  reactions,  to 
the  familiar  hydrodynamic  equations 
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(the  notations  hero  agree  with  those  of  L.  D.  Landau  [  1  ] )  w <•  must 
add  K  equations  of  continuity  for  each  chemical  component  K,  (l  = 

i 

-  1,  K)  participating  In  R  reactions,  H  (s  -  i,  R).  The 

following  I'orm  ol‘  these  equations  Is  used  below: 
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whore  c,  lr  a  portion  >f  the  mass,  1,  Is  dlffusi  .n  flow,  m,  Is  the 
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m  iocu.ar  weight  [g/molc]  1‘  the  component  K;,  and  then  q> .  Is  the 

velocity  .f  R  (mole/cm3  sec]  and  v.  Is  stoi ch lame  trie  coelTlcioi.t 
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V  v  no-dlmenslonal  steady  flows,  a  single  integration  -f  the 
extracted  equations  yield  (in  Landau's  notations  [l]) 
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where  C  with  the  subscripts  are  constants  of  Integrations.  We  shall 

.introduce  for  consideration  the  expression  for  thermal  flow  [1] 
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and  the  dimensionless  variables  . 
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Considering  the  thermal  capacities  constant,  we  represent  (5) 
and  (6)  in  the  dimensionless  form 
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change  the  pressure  1  r»  front  ..f  tne  name  fr-nt,  we  go 1 
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If  is  essential  that-  when  P  >  1,  the  pressure  In  !'r  nt  I'  the 
flame  !':•  :.s  It. crease;;  on  appr  inch  t.  >  tne  fr  a.t  rind  ■  •  a. verse  ly,  wner. 

P  •.  1,  the  pressure  decreases.  Note,  that  an  increase  Ir.  tne  pressure 
In  fr  a.t  f  the  flame  fruit  whet.  I’  >  1  an  be  a  s  .ur  e  <f  cl  Isiuri  at.  v . 
K  r  tne  state  It.  the  burn 'rtf,  save,  we  have  It.  ;  lie  sere  apprjxl- 
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•U  /<*,.  / 1,,—1  \ 

I  1 


and  !:.  the  first  approximati  t. 


These  results  are  obtained  provided  P  =  1,  and  p.  =  1.  for  e  )ee- 
tr  .ehc-mi  ■  a i  reactions  it  is  evident  that  <  0. 

Thus,  if  in  front  of  the  flame  front  we  have  ^  /  >  0  or 
i  ■'  0  for  P  >  1  and  P  <  1,  respect! ve.ly,  and  is  proportional  to 

Mi,  then  in  the  region  of  the  front  <  0  and  is  proportional  to 

Mi2  (q/°i2)  (<1  is  the  thermal  effect  of  the  reaction  fcal/fj);  here 
Mi2  «  1  and  q/cj  ~  10,  if,  for  example,  Q  ~  104  cal/mole,  p  ~  5^ 

ir 

g/mole,  and  cx  ~  JOO  m/sec. 

IiCt  us  now  consider  the  nonstationary  system  of  equations  out¬ 
side  of  the  flame  front.  It  Is  more  convenient  to  use  in  place  of 


Eq.  (3)  the  equivalent  form: 


du>  dp  (^ri 

p  7/7  —  at  3'i*  /ar-'l‘'"/ 


So 


Considering  an  before  the  thermal  capacity  of  the  medium  to  be 
constant,  let  us  examine  the  last  equation  In  a  zero  approximation 


(with  respect  to  M2) 


k, 


T  (V.  V)j  /■  div(xV/ , 


Combining  this  equal.  1  n  with  the-  equation  of  continuity  (1),  and 
by  using  the  equation  >  f  state,  where  p  =  const,  we  get 


Av  V— v,.,  (x  ,  r  )  VT 


(15) 


(it.  the  cylindrical  and  spherical  cases  Vi,.^  =  0). 

in  the  same  approx  I ma  t  J  - n  the  equation  for  determining  p  for 
centrally  symmetrical  flows  has  the  form 


■V,  f  r’W.,,  ,  1  "  1)  0 

I  i>r  \  c/,  t)r  (.  / 


(16) 


where  n  =  0,  1,  2,  respect! ve ly,  in  the  plane,  cylindrical,  and  spherical 
case.  Thus,  the  calculation  of  these  flows  in  the  zero  approximation 
reduces  to  a  solution  of  one  nonlinear  parabolic  equation,  for  which 
it  Is  still  necessary  to  formulate  specific  initial  and  boundary  con¬ 
ditions  . 

We  turn  now  to  the  calculation  for  changing  the  pressure,  limit¬ 
ing  ourselves  to  the  one-dimensional  case.  The  combination  of  Eqs . 

(l)  and(2)  yields  the  expression 


i)-;i  _  if1 

~  (l?‘ 


:,-[/>  +  prs-(  -I-  t)£J 


Hence 


p  ~  ( .T  1  +  ?)  iu  ~  1)1,1  +  \  \  Hx  +  const 


In  the  dimensionless  variables  already  used,  and  also  in  the  new 


SC  p  V 

J.  (It  ,  V  =*  —  ,  a  =  — 
M  P I i 3 


s/ 


ones 


Instead  of  (15)  and  (In)  we  wli.  have,  respect lvely 

i)0101 


Jv1"1  ,  .»  ,  ... 

-  =3  —  In  (|  .1.  v101) . 
i(t  ><X  i<X*  ' 


After  appropriate  substitutions  In  (17),  we  obtain 


1,111  T  (/*  —  1) 


‘7 


(18) 


it  can  be  shown  that  H  r  centrally  symmetric  flows  the  corre¬ 
sponding  equation  f  r  j_  1  rj  the  functions  p  and  v^  has  the  form 

dp  //,ll  I  .  x  \  ii  f  |  ,1  I  ! 

"’■ioK-  l~-f: 

A  comparison  <f  (18)  and  (1  )  with  consideration  of  the  equa- 
tions  defining  O'  '  in  b  -th  cases,  leads  us  to  the  conclusion  that 
when  the  Initial  temperature  distribution  coincided  with  stationary 
steady  solution  in  the  case  of  a  subsequent  Increase  in  temperature, 
unsteadiness  increases,  with  respect  to  the  modulus,  the  change  in 
pressure  in  front  of  the  flame  fr  >nt  (because  J^/d t  >  0).  The  sign 
of  pressure  change  is  usually  determined  by  the  sign  of  the  inequality 
P  £  1.  Unlike  the  stationary  case,  the  states  behind  the  flame  front 
here  can  be  described  not  only  by  identity  solutions. 

If  the  temperature  distribution  is  a  monotonically  increasing 
function  of  the  coordinate,  Its  curvature  behind  the  flame  front  is 
negative  and  the  sign  of  the  change  in  pressure  behind  the  front  is 
the  inverse  of  that  ahead  of  the  flame  front.  From  considerations  of 
continuity,  the  pressure  distribution  when  P  >  1  should  In  this  case 
have  a  maximum  at  the  leading  edge  of  the  flame  front  and  a  minimum 
at  the  trailing  edge;  for  P  <  1  the  pressure  distribution  will  be  a 
monotonically  decreasing  function  of  the  coordinate  (in  any  case 
exterior  to  the  flame  front).  The  density  distribution  as  evident 
from  the  equation  of  state,  when  p  =  const,  will  be  monotonically 
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decreeing  function  of  the  coordinate  (in  any  case  outside  the 
flame  front).  From  (lS)  we  see  that  here  the  velocity  in  a  zero 
approximation  should  pass  through  the  maximum  at  the  flame  front. 

If  the  temperature  distribution  has  a  maximum  at  the  flame  front, 
Its  curvature  outside  the  front  at  a  sufficient  distance  from  It  are 
positive,  and  the  changes  in  the  sign  of  the  pressure  in  front  of  and 
behind  the  front  are  identical.  When  P  >  1,  the  pressure  distribution 
should  have  a  maximum  at  the  leading  edge  of  the  flame  front,  and 
f  r  P  <  1  it  should  have  a  minimum  at  the  trailing  edge.  The  density 
distribution  should  have  a  minimum  at  the  flame  front,  and  the  velocity 
distribution  will  be  a  monotonlcally  Increasing  function  of  the  co¬ 
ordinate  (in  any  case  outside  the  flame  front). 

To  Illustrate  the  obtained  results,  Figs.  1-4  show  qualitative 
graphs  of  the  dimensiunless  characteristics  of  the  flow  (the  tempera¬ 
ture  0,  pressure  ?r,  density  v,  velocity  o)  as  a  function  of  the 
dimensionless  coordinate  X-  The  numerals  I,  II,  and  III  on  the  graphs 
correspond  to  the  regions  up  to  the  flame  front,  the  flame  front, 
and  behind  the  flame  front,  respectively.  Letters  A  and  B  refer  to 
the  two  possible  temperature  distributions:  the  montonically  increas¬ 
ing  distribution  (a)  and  the  distribution  with  the  maximum  at  the 
flame"  front  (b).  The  numbers  1,2  refer  to  the  cases  P  >  1  and  P  <  1. 

For  discussing  this  work,  we  are  Indebted  to  L.  N.  Prya tnltskiy, 
who  also  directed  the  author's  attention  to  the  possibility  of  increas¬ 
ing  the  pressure  at  the  leading  edge  of  a  flame  front. 
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